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Implementation of an optically active material on silicon has been a persistent technological challenge. For tandem 
photovoltaics using a Si bottom cell, as well as for other optoelectronic applications, there has been a longstanding 
need for optically active, wide band gap materials that can be integrated with Si. ZnSiP 2 is a stable, wide band gap 
(2.1 eV) material that is lattice matched with silicon and comprised of inexpensive elements. As we show in this paper, 
it is also a defect-tolerant material. Here, we report the first ZnSiP 2 photovoltaic device. We show that ZnSiP 2 has 
excellent photoresponse and high open circuit voltage of 1.3 V, as measured in a photoelectrochemical configuration. 
The high voltage and low band gap-voltage offset are on par with much more mature wide band gap III-V materials. 
Photoluminescence data combined with theoretical defect calculations illuminate the defect physics underlying this 
high voltage, showing that the intrinsic defects in ZnSiP 2 are shallow and the minority carrier lifetime is 7 ns. These 
favorable results encourage the development of ZnSiP 2 and related materials as photovoltaic absorber materials. 


Broader Context 

Of all the renewable energy technologies, solar photovoltaic electricity has one of the highest resource potentials; 
there is enough energy in the sunlight incident on the surface of the earth to meet the world’s energy demands 
many times over (~10,000:1). However, significant market penetration requires photovoltaics to be cost competitive 
with fossil fuels, even when unsubsidized. Currently, balance of system costs, rather than module costs, represent 
the majority of the total installed cost. Thus, increasing module efficiency is attractive as high efficiency cells can 
reduce installation size and therefore cost. Tandem photovoltaic architectures can provide a transformative boost 
in module efficiency over the single junction alternative due to reduced thermalization losses. Silicon photovoltaics 
is a well established (>90% market share), high efficiency, low cost technology that provides a crystalline template 
to grow top cells upon. However, the top cell material must satisfy strict requirements, including high efficiency and 
long reliability, or its presence will simply reduce the performance of the silicon bottom cell. The primary top cell 
materials considered to date include III-V materials, but the cost of these materials and their sensitivity to defects 
have proven challenging. In this work, ZnSiP 2 emerges as a wide band gap absorber that has the potential to meet 
the requirements needed for a top cell in tandem silicon-based photovoltaics. 
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Fig. 1 Theoretically determined band gaps versus lattice 
constants for some of the more earth abundant II-IV-V 2 
chalcopyrites. Also shown are III-V materials along with Si 
and Ge from group IV. The gray vertical bars overlaying Si 
and Ge highlight materials with similar lattice constant, 
within ± 1%. 


1 Introduction 

Optoelectronic materials are typically semiconductors 
that are covalently bonded with tetrahedral coordina¬ 
tion (i.e. group IV, II-VI and III-V materials). While 
the ternary II-VI analogs such as I-III-VI 2 (e.g. CuInSe 2 ) 
have been extensively investigatedit is surprising that 
III-V analogs such as II-IV-V 2 (e.g. ZnSiP 2 ) have re¬ 
ceived less attention for optoelectronic devices.®^ Fig¬ 
ure [T] shows that II-IV-V 2 chalcopyrite compounds span 
a broad range of band gaps, from 0.4 eV (CdSnAs 2 ) to 
2.5 eV (MgSnP 2 ), and have lattice constants that are 
compatible with many unary and binary semiconductors. 
Many of the II-IV-V 2 compounds are formed from com¬ 
paratively inexpensive and non-toxic elements (e.g. Zn, 


Mg, Si, Sn, N, P), rendering them particularly attractive 
for applications requiring large scale deployment such as 
photovoltaics (PV).®®Wide band gap II-IV-V 2 materi¬ 
als have seen virtually no investigation, despite significant 
opportunities for tandem PV, light emitting diode, pho¬ 
tonic circuit, and laser applications.^®^ Within tandem 
PV, it has been challenging to find wide band gap materi¬ 
als suitable for tandems, and particularly those compat¬ 
ible with Si. 17 ilSi Thus, the II-IV-V 2 compounds ZnSiP 2 
and ZnGeP 2 are of great interest as epitaxial top cell ma¬ 
terials on a silicon bottom cell. Implementing these mate¬ 
rials as inexpensive, earth-abundant top cells on Si lever¬ 
ages the dominance of Si PV (>90% market share) .®The 
work presented herein is focused on PV relevant charac¬ 
terization of ZnSiP 2 - 

The fundamental properties of ZnSiP 2 have been stud¬ 
ied since the late 1950’s using crystals that have been 
grown in a flux (typically Zn or Sn) or by halogen as¬ 
sisted vapor transport.® 11 12 20 37 Studies of these crys¬ 
tals reveal that ZnSiP 2 has a very small lattice mismatch 
with Si of 0.5% (Fig. [I]), has a band gap of ~2.1 eV, 
forms with minimal atomic disorder, and is structurally 
stable at temperatures up to 800 °Ci 1 4 1 -® ® Doping of 
ZnSiP 2 has yielded n-type (Se, Te, In, or Ga) and p-type 

(Cu) crystals.™™!! 

Some characterization has been done which is specif¬ 
ically relevant to the applications of tandem PV cells 
with silicon. Several authors have proposed and dis- 
cu ssed the prospects of ZnSiP 2 heterojunctions with 
Si. l 12 l 26 l 4 iJ42j Growth of Si/ZnSiP 2 interfaces has been 
demonstrated through heteroepitaxial crystallization of 
Si on ZnSiP 2 substrates,® growth of polycrystalline 
ZnSiP 2 on Si®, and epitaxial ZnSiP 2 on Si by vapor- 
liquid-solid growth.® While photoconductivity has been 
dem onstr ated, no PV devices have been realized to 
date.® 3 ® A ZnSiP 2 /Si device may be expected to have 
good light transmission through the top cell and into the 
bottom cell, because ZnSiP 2 has little parasitic below- 
band-gap absorption and good index of refraction match¬ 
ing with Si (reflection at the Si/ZnSiP 2 interface would 
be less than 1%).® Considering that it is a stable, bipolar 
dopable, wide band gap semiconductor that is lattice and 
index of refraction matched with Si, ZnSiP 2 is a promis¬ 
ing material for implementation as a top cell on Si PV. 

In this work, we have employed a variety of theoreti¬ 
cal and experimental techniques to address some of the 
challenges and unknowns regarding ZnSiP 2 as a top cell 
material on Si PV. Through a combination of photolumi¬ 
nescence measurements and first principles calculations, 
we show that the intrinsic material contains both donor 
and acceptor defects, and the associated energy levels are 
shallow (~0.1 eV from the associated band edge), result- 
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parameters were allowed to refine, the occupancies of all 
sites remained stable at 1.00(2) (see Supplementary Ta¬ 
ble S2 1 , indicating that all sites are fully occupied. A 
structural model with Zn and Si occupying both 4a and 
4b sites was also employed. This model led to less than 
1.3% occupancy of Si on the 4a site and no mixing of 
Zn onto the 4b site. Statistics for this model are com¬ 
parable to statistics for the model with no site mixing. 
In addition, all atomic displacement parameters led to 
nearly spherical thermal ellipsoids. Consistent with pre¬ 
vious results, and in contra st to other ternary and quar¬ 
ternary PV materials, I 49 ! 52 * ® ZnSiP 2 prefers a highly or¬ 
dered atomic structure that is very stable. ® 38 l 39 l 


Fig. 2 (a) Photograph of ZnSiP^ single crystals grown in 
Zn flux, (b) Crystal structure of II-IV-V2 chalcopyrites 
(I 42 d) emphasizing two interlaced networks of corner sharing 
tetrahedra. In one network, the tetrahedra (red) have a 
Group IV atom at the center and Group V atoms on the 
corners. The second network differs from the first, in that 
the tetrahedra (blue) are slightly distorted and they have a 
Group II atom at the center. 

ing in defect-tolerant behavior. We identify the origin of 
these defects as antisite defects. We show photoelectro¬ 
chemical (PEC) measurements of ZnSiP 2 which demon¬ 
strate excellent photoresponse and high open circuit volt¬ 
age (V oc ) of 1.3 V. Despite this being the first report 
of a ZnSiP 2 photovoltaic device, the observed voltage is 
higher than any other material lattice matched to Si, to 
our knowledge.®®! These findings establish ZnSiP 2 as a 
potential monolithic top cell material on Si. 

2 Results and Discussion 

2.1 Synthesis and Structure 

Flux synthesis yielded ZnSiP 2 crystals which are translu¬ 
cent, dark red in color, and have dimensions up to 2 
mm thick, 5 mm wide, and 20 mm long (see Fig. 1(a)). 
Single crystal X-ray diffraction (SCXRD) data confirms 
that ZnSiP 2 crystallizes in the I42d space group (No. 
122) with unit cell dimensions of a = 5.3986(2) and 
c = 10.4502(6) A (see Fig.[2](b)). The SCXRD structural 
analysis results are given in Supplementary Tables |S3f 
|S5| and agree well with previously reported structural 
data.® 

Some chalcopyrites are known to display atomic dis¬ 
order resulting in variation of their optical proper¬ 
ties.^ 23 ®^ Refinements of the SCXRD data indicate no 
significant disorder of Zn and Si between the 4a (Zn) and 
4b (Si) sites. When occupancies and atomic displacement 


2.2 Photoluminescence Characterization 

A combination of photoluminescence (PL) measurements 
and first-principles calculations (Section 2.4) of intrin¬ 
sic point defects (vacancies, antisite defects) were used 
to understand the recombination mechanisms in ZnSiP 2 - 
Photoluminescence measurements were performed across 
a range of temperature (5-300 K) and excitation power 
(0.1-15 mW). A large PL peak at 1.8 eV (Fig. [t] (a)) 
was observed at all measurement temperatures. Con¬ 
sistent with prior work, a second small sharp peak at 
2.04 eV was observed at temperatures below 15 K and 
is attributed to excitons.® The time dependence of the 
1.8 eV peak was measured to further probe the minority 
carrier recombination rates. In concert, these measure¬ 
ments suggest donor-acceptor pair (DAP) recombination 
via shallow defects. 

Figure |3] (b) shows the time dependence of the PL in¬ 
tensity associated with the 1.8 eV, at room temperature. 
An approximation of the carrier lifetime can be achieved 
using a biexponential model of the form 

I{t) = Ae- t/T ‘ + Be~ t / Tl , ( 1 ) 

where A , B 1 r s , and 77 are fit parameters (orange curve 
in Figure [3] (b)). This model assumes there are two inde¬ 
pendent, dominant decay mechanisms, and that each be¬ 
haves exponentially. The parameters r s and 77 are short 
and long lifetimes, respectively, and were found to be 
t s = 0.9 ns and 17 = 7.1 ns. This lifetime is compara¬ 
ble with those of many well developed thin film photo¬ 
voltaic m ateria ls such as CdTe (~20 ns) and Cu 2 ZnSnSe 4 
(~18 ns).®® From prior work, the hole mobility, /j, p , has 
been measured as high as 25 cm 2 /Vs.®® For these n- 
type crystals, a hole diffusion length (at the surface) of 
680 nm is calculated from /r p and 77 . 

If DAP transitions are the primary PL mechanism, the 
recombination kinetics depend on the wavefunction over¬ 
lap between the spatially separated electron, trapped on 
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Fig. 3 (a) Photoluminescence spectrum of ZnSiP 2 (red) and absorption edge (grey). The inset energy level diagram 
illustrates the most likely mechanism for the observed luminescence: donor/acceptor pair (DAP) transitions, (b) Time 
dependence of the 1.8 eV peak is evidence of DAP transitions and shows good carrier lifetime, (c) The power dependence of 
the 1.8 eV peak position shows a blue shift characteristic of DAP transitions, (d) The power dependence of the integrated 
intensity of the 1.8 eV peak; a fit to the data with I(P) oc P k further supports DAP transitions, (e) The temperature 
dependence of the integrated intensity of the 1.8 eV peak fits Eqn. (|5j), giving Ed = 58 meV. 


a donor site, and hole, trapped on an acceptor site. 57 ® 
As such, the PL recombination is more complex than the 
two-process model assumed aboveJ 57 ! The recombination 
rate for DAP transitions as a function of the distance, r, 
between donor and acceptor is assumed to be of the form 

W{r) = Wq exp [—2r/a \, (2) 

where Wq is the rate as r —>• 0 and a is the larger of the 
donor or acceptor effective Bohr radii. 57 The PL intensity 
as a function of time, t , is proportional to 


pOO 

I(t)ccexp 4 tt N ( exp[—W{r)t} — \)r 2 dr 

V Jo 

POO 

x 47 t N W(r) exp[—W(r)t\ r 2 dr, (3) 
Jo 


where N is the concentration of the majority defect (pos¬ 
sibly donors in this case because PEC rectification in¬ 
dicates n-type conductivity).® While Wq, a, and N 
are free parameters, the parameters a and N are cou¬ 
pled into the non-dimensional constant £ = a 3 N, and 


Eqn. § and were fit to the data in non-dimensional 
form (see Eqn. ( |S1[ ) and (S2|). The black curve in Fig¬ 
ure [3] (b) shows that the resulting fit is significantly bet¬ 
ter than the biexponential fit (orange), despite fewer free 
parameters. Because the donor and acceptor concentra¬ 
tions are similar (as suggested by our theoretical analy¬ 
sis in section [2.3[ ), the donor and acceptor contributions 
to £ = (2.3 ± 0.3) x 10~ 2 cannot be decoupled. Wq 
was found to be 8 x 10 9 s _1 , slightly higher than GaAs 
(3 x 10 8 s" 1 ).® The fit of Eqn. ^ to the time depen¬ 
dent data (Fig. [3] (b)) demonstrates that DAP transitions 
are likely the dominant recombination mechanism. This 
conclusion is in contrast to previous analysis of time de¬ 
pendent cathodoluminescence done at 10 K rather than 
room temperature, which determined the recombination 
to be due to bound excitons with phonon sidebands.® 
This difference is likely due to different recombination 
mechanisms dominating at different temperatures. 


The power dependence of the PL gives additional evi¬ 
dence of DAP transitions. With increased injection level, 
the increased density of excited charge carriers leads to 
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a stronger Coulomb interaction between ionized DAP’s. 
This interaction leads to an increase in the emitted pho¬ 
ton energy following: 


hu(r) = Eg - E a - E d H-( cgs ), 


(4) 


where E g is the band gap, Ea is the ionization energy of 
the acceptor, Ed is the ionization energy of the donor, e 
is the electron charge, and r is the distance between the 
donor and acceptor responsible for the transition. Fig¬ 
ure [3] (c) shows the peak position shifting to higher en¬ 
ergy as the injection level is increased, characteristic of 
DAP transitions. 

Finally, the power dependence of the peak intensity 
(Fig. [ 3 ] (d)) can be considered within the DAP hypothesis. 
An exponential model of the form /(P) oc P fc , where P is 
the excitation power and k is the exponent, was fit to the 
experimental data. For the nine distinct temperatures at 
which power dependent data was collected, the average k 
was 0.72 with a standard deviation of the average of 0.12. 
With k < 1, the power dependence of the peak intensity is 
characteristic of either free to bound, or DAP transitions 
as being the primary source of radiative recombination. 61 

In addition to recombination mechanisms, PL data 
analysis can provide the donor activation energy (the ma¬ 
terial is n-type, as shown in Sec. 2.4). Fig. [ 3 ] (e) shows 
that the integrated intensity was found to increase signif¬ 
icantly with decreasing temperature down to 25 K. Below 
25 K, the PL intensity decreased and eventually appeared 
to saturate. The temperature region above 25 K was fit 
with 


J ( T ) ^ 1 + ae (E D /k B T)’ ( 5 ) 

where a is a process rate parameter, kp is Boltzmann’s 
constant, and T is the temperature.!- 62, Six temperature- 
dependent data sets were collected at different excitation 
powers, and the activation energy of the 1.8 eV peak was 
found to be Ed = 58 ± 16 meV; this value is shallow, and 
in agreement with previous reports] 60 * 63 ^ 

Both the time and power dependence of the PL in¬ 
tensity and the power dependence of the peak position 
indicate DAP recombination. There was no evidence in 
the PL of deep level defects that would be detrimental 
to device operation, and the temperature dependence of 
the PL intensity predicts the donor level to be shallow, 
corresponding with the high V oc found in the PEC char¬ 
acterization. 

2.3 Identification of Potential Intrinsic Defects 

To further understand charge carrier recombination in 
ZnSiP 2 , first-principles density functional theory (DFT) 


calculations were employed to identify possible point de¬ 
fects that act as donors and acceptors in ZnSiP 2 - By 
calculating the enthalpies of formation of known Zn-Si-P 
phases, a stability map in chemical potential space was 
constructed (Fig. [ 4 ] (a)).^ The enthalpy of formation 
of intrinsic defects within the ZnSiP 2 region of stabil¬ 
ity of this chemical potential space were also calculated 
(Fig. m (b)), along with the associated density of states 
(DOS) of relevant defects (Fig. [4] (c)). The formation 
enthalpies (AHD, q ) of 7 different point defect types (3 
vacancies, 4 antisite defects) were calculated in 7 differ¬ 
ent charge states (3- through 3+; 4- and 4+ were addi¬ 
tionally examined for Vsi) as a function of Zn, Si and P 
chemical potentials, and Fermi level (Ep). Although pre¬ 
vious electron paramagnetic resonance experiments have 
shown evidence of Zn and P vacancies, as well as Si on 
Zn antisite defects,! OHM] these are the first calculations 
identifying likely point defects as a function of chemical 
potential for ZnSiP 2 - 

The Zn-Si-P phase diagram calculated at standard con¬ 
ditions using Fitted Elemental-Phase Reference Energies 
(FERE) is shown in Fig. [ 4 ] (a) ZnSiP 2 occupies a large 
region of this chemical potential phase space, suggesting a 
large window for synthesis conditions. Competing phases 
include Zn, Si and P elements, and the binary compounds 
SiP, Zn 3 P 2) and ZnP 2 . Within the single phase region, 
where ZnSiP 2 is stable, there are three distinct regimes 
(labeled Rl, R2, and R3) with different lowest <5 L/]d g de¬ 
fect pairs. The dominant DAP’s that set the Fermi level 
are Si|J and Sip” in Region 1, Pg(“ and Zn|r in Region 

2, and SiijA and Zn|[~ in Region 3. The Fermi levels {Ep) 
are indicated by dashed vertical lines in Fig. 0(b)- 

The ZnSiP 2 crystals were grown under Zn-rich, Si- 
poor, and P-poor conditions. These growth conditions 
most likely lie near the line connecting points 1, 2, and 

3, shown in Fig. [4] (a). For this reason, the defect forma¬ 
tion enthalpies and their concentrations were examined 
at these points which lie in the corresponding regions Rl, 
R2, and R3, in Fig. [4] (a). For the point defects consid¬ 
ered, AHu,q as a function of Ep is shown in Fig. 0(b) 
wherein the formation enthalpy of the lowest-energy DAP 
at each of the three points is shown in bold. A first order 
approximation of the defect concentrations can be deter¬ 
mined from 


[C] = [5] e" 


Af* D ,q 
k B T 


( 6 ) 


where [S'] is the crystallographic site concentration 
([Zn] = [Si] = [P] /2 ss 1.3 x 10 22 cm -3 ), kp is the Boltz¬ 
mann constant, T is the temperature at which the defects 
form, and ALAyq is evaluated at the chemical potentials 
associated with a particular point in Fig. [4] (a) and the 
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Fig. 4 (a) Zn-Si-P phase diagram in chemical potential phase space. The flux growth conditions are Zn rich and are therefore 
expected to lie in the region between point 1 and point 3. Also shown is a heat map representing the Fermi level (Ef), which 
ranges from near mid-gap at point 1 to ~1.7 eV above the valence band maximum. Rl, R2 and R3 are regions in the chemical 
potential phase space where different donor-acceptor pairs are dominant, (b) Defect formation enthalpies (AH_d, 9 ) of 7 
different point defects (vacancies and antisites) in all possible charge states (q ranging from 4- to 4+) at points 1, 2, and 3 on 
the phase diagram (a). The Fermi level is shown as a vertical dashed line along with mid-gap which is shown as a vertical 
dotted line; the relevant defects which set the Fermi level are labeled in bold. For additional details on all charge states of 

(c) The density of states (DOS) of ZnSiP 2 , without defects, and ZnSiP 2 with four key native 
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Table 1 Concentrations of dominant point defects in 
ZnSiP 2 at 500 °C. The locations of points 1, 2, and 3 are 
identified in Fig. [4] (a). [D] and [A] are the donor and 
acceptor concentrations in cm -3 . 


Point 

DAP 

[D] 

[A] 

1 

Sizi/Sip 

7.5 x 10 18 

1.5 x 10 19 

2 


6.1 x 10 13 

6.1 x 10 13 

3 

Ps+/Zn|r 

4.3 x 10 15 

4.3 x 10 15 


corresponding Fermi level. The defect concentrations 
were calculated at T = 500 °C, resulting in a lower limit 
to what may be observed experimentally; this is because 
the crystal growth occurs as the temperature of the Zn 
flux slowly cools (-2 °C/h) from 1000 °C to 500 °C, 
after which the flux is rapidly cooled to room tempera¬ 
ture. The calculated DAP’s and their associated concen¬ 
trations are given in Table [l] 

The density of states (DOS) of defect-free ZnSiP 2 and 
ZnSiP 2 containing one of the four prominent defects are 
shown in Fig[4](c). The shallow donor and acceptor states 
appear as shoulders near the conduction band minimum 
(CBM) and the valence band maximum (VBM), respec¬ 
tively. There is no evidence of mid-gap states in any 
case. Previous studies have observed a deep acceptor 
level, ^0.7 eV above the VBM, which appears to have 
been the result of extrinsic impurity (Cu) doping.^ 1 * 24 * 67 ! 

PL experiments and theoretical defect calculations are 
consistent with DAP antisite defects as the prevalent re¬ 
combination mechanism in ZnSiP 2 . These DAP defects 
work together to set the Fermi level such that ZnSiP 2 is 
expected to be intrinsic, or moderately n-type, across a 
wide range of synthesis conditions. Analysis of the tem¬ 
perature dependence of the PL intensity found the donor 
level to be shallow: ~58 meV below the CBM. Corre¬ 
spondingly, the dominant defects identified in the cal¬ 
culations were found to occupy shallow levels, based on 
their associated DOS. These results indicate that, over 
the full chemical potential region of stability, ZnSiP 2 PV 
devices are defect-tolerant and should have a high V oc , in 
agreement with the PEC results from Section [274} 

2.4 Photovoltaic Characterization 

Regenerative photoelectrochemistry (PEC) was used to 
characterize the photovoltaic properties of ZnSiP 2 . This 
technique allows the measurement of the current den¬ 
sity vs. potential ( J-E ) of irregularly shaped bulk sin- 
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Fig. 5 (a) Schematic of the photoelectrochemical cell with 
a ZnSiP 2 working electrode and Si photodiode for 
calibrating light intensity, (b) Photoelectrochemical J — E 
data for ZnSiP 2 electrodes under 1 Sun illumination in 
contact with 50/350 mM Fe(CN)g 3 ^ ■* (blue), and 5/5 mM 
Fe(CN)g 3 1 (light=green, dark=black), (c) Spectral 
response of the same (blue in (b)) ZnSiP 2 electrode in 
50/350 mM solution, along with the percent absorption of 
3 mm of solution and of the 0.5 mm thick ZnSiP 2 electrode. 
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gle crystals without the need to create a solid state p-n 
junction. A ZnSiP 2 crystal is formed into an electrode 
and submersed in an electrolyte solution containing a 
one-electron, reversible, outer-sphere redox couple, cre¬ 
ating a semiconductor-liquid junction (Fig [5] (a)). This 
method has proven to be a successful initial characteri¬ 
zation technique for many semiconductors.®®^ Within 
this photoelectrochemical configuration, initial measure¬ 
ments as high as 0.9% photoconversion efficiency were 
obtained from ZnSiP 2 single crystals in contact with the 
Fe(CN)g 3 ^ (aq.) redox couple. 

Photoelectrochemical electrodes were constructed from 
several single crystals, and representative J-E data 
is shown in Fig. [5] (b). The representative J- 
E data was collected for electrodes in contact with 
50 mM K 3 Fe(CN) 6 /350 mM K 4 Fe(CN) 6 (aq) and 5 mM 
K 3 Fe(CN) 6 /5 mM K 4 Fe(CN) 6 (aq).® The photoelectro¬ 
chemical rectification with the Fe(CN)g 3 ^ (aq) elec¬ 
trolyte indicates that the ZnSiP 2 crystals are n-type. 

In the 50/350 mM Fe(CN)g - ^ 4- electrolyte, an E oc 
of 1.08 V, a short circuit current density (J sc ) of 
0.73 mA/cm 2 , and a fill-factor (FF) of 0.53 were mea¬ 
sured. However, at these concentrations, the redox cou¬ 
ple effectively absorbs all incident AM 1.5 light at ener¬ 
gies above 2.6 eV, which is a significant fraction of the 
illumination above ZnSiP 2 ’s 2.1eV optical band gap (see 
Fig E (c)). With a reduced concentration of 5/5 mM 
Fe(CN)g - ^ 4- , better E oc and J sc were obtained due to 
reduced absorption losses. At this lower concentration 
of the redox couple the champion electrode produced an 
E oc of 1.30 V, a J sc of 1.72 mA/cm 2 , a FF of 0.41, and 
an efficiency of 0.9%. This unoptimized cell configura¬ 
tion is limited by solution resistance, contact resistance, 
and concentration overpotential, each of which could be 
significant in these preliminary electrochemical measure¬ 
ments. The J-E data has not been corrected for mass 
transport limitations or concentration overpotential, and 
thus the low FF does not reflect the intrinsic performance 
of the semiconductor.® 

High E oc (Kc in a solid-state device) is indicative of 
good material quality, good charge carrier lifetime, and a 
lack of deep level defects. An ideal homojunction ZnSiP 2 
PV device with optimized doping of the base and emitter 
would be expected to have a maximum V oc of ~ 1.8 V jTfillll 

The Fe(CN)g 3 ^ 4 redox couple is not optimized for the 
band edge positions of ZnSiP 2 , nor has the effect of dop¬ 
ing of the ZnSiP 2 crystals been investigated. Therefore, 
the observed 1.30 V E oc of ZnSiP 2 under 1 sun condi¬ 
tions is promising because increasing the V oc can be one 
of the most significant challenges in the development of a 
new material.® The high V oc and low band gap-voltage 


offset observed here are on par with much more mature 
wide band gap III-V materials in solid-state devices E§HH1 
but in a material that can be easily integrated with Si. 
For materials lattice matched to Si, much lower voltages 
are typically observed, even in more optimized solid-state 
devices. For example, GaNPAs devices on Si have a re¬ 
ported V oc of 1.1 V, as well as poor fill factor and cur¬ 
rent.® A leading contender for earth-abundant tandem 
PV on Si is CZTS, but these devices on Si have a reported 
V oc of only 0.12 V,® and optimized devices not grown 
on Si have achieved only 0.66 V.® Thus, it is extremely 
promising that in this first report of photovoltaic perfor¬ 
mance of ZnSiP 2 , we have observed such a high voltage. 


There are several potential explanations for the ob¬ 
served low current density of the electrodes: diffusion 
limitations due to crystal thickness, weak light absorp¬ 
tion, resistive losses, and parasitic light absorption by the 
solution. To better understand the photocurrent gener¬ 
ation in these devices, the electrochemical spectral re¬ 
sponse was investigated. Figure [5] (c) shows the external 
quantum efficiency (EQE) of a typical ZnSiP 2 electrode 
in 50/350mM Fe(CN)g - ^ 4- (aq) along with the percent 
absorption by the solution (UV-Vis), and the percent ab¬ 
sorption by the ZnSiP 2 (photothermal deflection spec¬ 
troscopy®). The J-E performance of the same electrode 
under AM1.5G illumination is shown (blue) in Fig [5] (b). 
From the spectral response, it is apparent that the short 
wavelength spectral region (> 2.6 eV) is strongly ab¬ 
sorbed by the electrolyte so only a small band of the 
AMI.5 spectrum (^2.1-2.6 eV) is absorbed by the crys¬ 
tal. The long wavelength edge of the spectral response 
closely resembles the ZnSiP 2 absorbance spectrum, indi¬ 
cating that photocurrent is limited by light absorption, as 
expected due to the symmetry forbidden nature of transi¬ 
tions at the band gap J 29 * 38 ® This arises because the elec¬ 
tric dipole transition between the valence band maximum 
(VBM) and conduction band minimum (CBM) is forbid¬ 
den by a selection rule.® The crystals were also thick 
(0.5-2 mm) relative to the expected absorption depth of 
the photons, so it is likely that many photocarriers re¬ 
combined in the bulk before they could be collected. The 
total J sc , as calculated from integrating the spectral re¬ 
sponse, was 0.35 mA/cm 2 , 27% lower than that measured 
in the J-E characteristic (Fig. [5] (b)). While the spec¬ 
tral response and J-E characteristic were collected from 
the same ZnSiP 2 electrode, the experimental setups were 
different in each case. These differences are primarily 
responsible for the discrepancy in J sc between the two 
different measurements. 
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3 Conclusions 

A combination of photoeletrochemical (PEC) charac¬ 
terization, photoluminescence (PL) measurements, and 
theoretical intrinsic defect calculations have shown that 
ZnSiP 2 has the potential to be a high quality photo¬ 
voltaic top cell material. PEC characterization revealed 
an open circuit voltage of 1.3 eV, indicating a lack of 
deep level defects. PL confirmed the presence of shal¬ 
low donor-acceptor pairs as well as the lack of deep level 
defects. Theoretical defect calculations established that 
the dominant defects in ZnSiP 2 are shallow antisite de¬ 
fects throughout the stability window in Zn-Si-P chemi¬ 
cal potential space. ZnSiP 2 is additionally predicted to 
be intrinsic, or moderately n-type throughout this entire 
region in chemical potential space. Thus, this material 
is intrinsically very stable and defect-tolerant. In this 
first report of the solar energy conversion properties of 
ZnSiP 2 , we have shown that it has excellent potential as 
a top cell on Si due to its abundance, defect tolerance, 
high voltage potential, and wide phase stability window. 

4 Methods 

4.1 Experimental 

Large single crystals of ZnSiP 2 (needle shaped or some¬ 
times hollow tubes, up to 2 by 5 by 20 mm) were syn¬ 
thesized using a zinc self-flux.® Zn (Alfa Aesar, #10759, 
99.999% metals basis), Si (Alfa Aesar, #43006, 99.9999% 
metals basis) and P (Alfa Aesar, #10670, 99.999% met¬ 
als basis) were mixed in 20:1:2.5 molar ratio and sealed 
in a 19 mm inner diameter quartz ampule under static 
vacuum better than le-6 mbar. To facilitate dissolu¬ 
tion in the flux, the silicon was ball milled (Spex 8000 
Mixer/Mill) in silicon carbide vials to a fine powder be¬ 
fore use. The sealed ampule was placed vertically in a 
tube furnace and heated to 1000 °C and soaked for 30 
hrs to ensure complete dissolution and homogenization. 
After the soak the furnace was cooled at 2 °C/hr to 500 
°C when it was removed. The excess molten flux was 
immediately decanted by tilting the hot ampoule before 
it cooled to room temperature. Residual Zn flux was 
removed from the crystals by etching with hydrochloric 
acid (37% by volume) mixed with deionized water in 1:3 
ratio by volume. Residual Si precipitates were removed 
from the crystal surfaces by etching in a dilute mixture of 
hydrofluoric (49% by volume) and nitric (70% by volume) 
acids in deionized water in 10:3:87 volume ratio. Other 
precipitates or inclusions have not been observed. The 
crystals did not receive further treatment before mea¬ 
surements were performed except when otherwise noted. 


A fragment (approximately 0.01x0.01x0.03 mm) cut 
from a larger single crystal was used to collect single 
crystal X-ray diffraction (SCXRD) data. This crystal 
fragment was mounted onto the goniometer of a Bruker 
KappaCCD diffractometer equipped with Mo K a (A = 
0.71073 A) radiation. Data collection and structure solu¬ 
tion were performed with SHELX-97. Structural refine¬ 
ments and extinction corrections were performed using 
SHELXL suite. 

The photoelectrochemical (PEC) response of ZnSiP 2 
was tested using the Fe(CN)g 3 ^ redox couple (potas¬ 
sium ferrocyanide and potassium ferricyanide in H 2 O 
with no supporting electrolyte).^® ZnSiP 2 electrodes 
were constructed by evaporating metal back contacts on 
one flat side of the crystals. The back contacts were 
700 nm of In, and they were capacitively blasted to 
achieve Ohmic behavior. A copper lead wire was affixed 
to the back contact with silver paint and the electrode 
was sealed at the end of a glass tube with epoxy (Loctite, 
Hysol 1C). An AM1.5G light source (ABET Technolo¬ 
gies, model no. 10500) was used to illuminate the front 
of the crystal and the light intensity was set by placing 
a calibrated photodiode (Thor Labs part no. FDS100) 
electrode in the same position as the crystal electrode, 
without electrolyte present, resulting in very repeatable 
data. Current-voltage data was collected using a Bio- 
Logic SP-240 potentiostat with Pt foil as the counter 
electrode and a saturated calomel reference electrode. 
Potentials were later adjusted to reference the solution 
potential, E{A/A~) : which was measured with a Pt wire 
in the electrolyte. 

Non-aqueous PEC data was taken using the same elec¬ 
trodes, which were briefly etched in dilute HC1, and then 
inserted into a 3-neck PEC cell on a schlenk line. The 
electrolyte consisted of 10 mM ferrocene (Fc, Sigma- 
Aldrich), and 0.05 mM ferrocenium dissolved in acetoni¬ 
trile with 0.1 M LiC 104 as a supporting electrolyte.® 
The same ABET model of light source was used for il¬ 
lumination and the photon flux was calibrated using a 
photodiode using a similar method to the aqueous tests. 

Photoluminescence spectra were collected with a 442 
nm HeCd laser with a 250 /iin spot size. Spectra were 
collected at temperatures of 5, 10, 15, 20, 25, 30, 50, 75, 
100, and 300 K. At each temperature, the laser power 
was varied and spectra were collected at 0.1, 0.5, 1.0, 5.0, 
10.0, and 15.0 mW (laser power at crystal surface). 

A Yb:KGW laser and an optical parametric amplifier 
(pulse length 0.3 ps, repetition rate 1.1 Mhz, wavelength 
420 nm) were used for single photon excitation photolu- 
minescence lifetime measurements. The excitation beam 
diameter was 0.2 mm and average laser power was 2 mW. 
Photoluminescence was measured with a 10 nm bandpass 






filter at 650 nm, a Si avalanche photodiode detector, and 
time-correlated single photon counting. Measurements 
were done at room temperature on two different, oppos¬ 
ing sides of the same crystal. The results were slightly 
different from one side to the other and measurements of 
different spots on the same side yielded very consistent 
data. 

4.2 Computational 

To study the intrinsic point defect physics in ZnSiP2, we 
performed density functional theory (DFT) calculations 
using the VASP codeJs3, in the generalized gradient ap¬ 
proximation (GGA) of Perdew-Burke-Ernzerhof,^ with 
the projector augmented wave formalism^. The total en¬ 
ergies were calculated for supercells containing 64 atoms 
(2x2x2 unit cells) using a plane-wave cutoff of 400 eV 
and a 4x4x4 Monkhorst Pack k-point grid. The defect 
supercells were relaxed following the numerical approach 
described in Ref. El for application of the electrostatic 
Hubbard U. Defect supercell total energies were calcu¬ 
lated for D = V Zn , V Si , V P , Zngj, Si Zn , Si P , Psi in charge 
states q = -3, -2, -1, 0, 1, 2, 3 for all defects and also q = 
-4, 4 for Vgi. Using the calculated total energies, the 
defect formation enthalpies (A Hp, q ) are calculated as 

A Ho, q = (E DtQ - E h ) + A 4 a + qE f (7) 

where the first term on right hand side is the energy dif¬ 
ference between the supercell with defect ( D ) in charge 
state q and the defect-free host (Eh)- The second and 
the third terms describe the chemical potential of the 
atomic and electronic reservoirs, respectively. The chem¬ 
ical potential /x a (a = Zn, Si, P) is given, relative to the 
elemental phase, such that + A/i q , where /i° is 

the elemental phase chemical potential, calculated using 
the Elemental Phase Reference Energies (FERE)® and 
A/i Q is the deviation in chemical potential from the ele¬ 
mental reference state. The Fermi energy (Ep) is given 
with respect to the valence band maximum (VBM) i.e. 
Ep = E v + AEp ; the value of Ep can vary from zero to 
the band gap. 

The underestimation of the band gap within the DFT- 
GGA method was “corrected” by applying band edge 
shifts, which are based on the GW quasi-particle energy 
calculations, as described in Ref. (SO! and [511 The band 


edge shifts for ZnSiP 2 were calculated to be A E V bm = 
-0.907 eV and A Ecbm = -0.084 eV. The band gap, af¬ 
ter application of the GW-calculated shifts, is 2.17 eV, in 
good agreement with experimental measurements.^ 2 2HE2 
While the general methodology for calculation of defect 
formation enthalpy from supercell calculations is well es¬ 
tablished, it must be noted that several corrections need 
to be applied to obtain a more accurate defect formation 
enthalpy. Following the methodology described in Ref. 
l93l we have applied the following corrections to A Hp> q -. 
(1) image charge correction for charged defects, (2) poten¬ 
tial alignment corrections for charged defects, (3) band 
filling corrections for shallow donors/acceptors, and (4) 
band gap corrections for shallow donors/acceptors. Some 
of these corrections (3, 4) depend on whether the de¬ 
fect is a shallow or deep level; we performed extensive 
tests, including DFT calculations with hybrid function¬ 
als (HSE06)®to ascertain the nature of the defect levels. 
To calculate the Fermi level, as a first approximation, we 
used the crossing point between the donor and acceptor 
with the lowest formation enthalpies on a A Hp> q vs. Ei? 
plot shown in Fig [4] (b). The defect concentrations [C] 
were calculated by multiplying the total number of lat¬ 
tice sites [5] that can accommodate a certain defect, by 
the Boltzmann factor (Eqn. [6|. 
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Electronic Supplementary Information 
Summary of II-IV-V 2 properties 


Table SI Experimental electronic & optical properties of II-IV-V 2 compounds. Forbidden transitions are shaded in red. All 
band gaps are for chalcopyrite phases. Experimentally determined chalcopyrite lattice parameters (a & c) are in A, band gaps 
(E“ p from experiment and E c from our GW calculations) in eV, mobilities ( p ,) in cm 2 /Vs and free carrier concentrations 
(n & p) in cm -3 . 


ID 

Material 

a 

c 

c/a 

E exp ( E caic) 

He 

n 

Hh 

V 

References 

0 

Si 

5.43 



1.12 

c 1500 

10 14 -10^ 

< 500 

10 14 -10 19 

"1551 ‘ 

1 

ZnSiP 2 

5.40 

10.44 

1.93 

2.0-2.3 (2.09) 

50-1,000 

10 13 -10 18 

1-25 

<10 T7 

TTTJ^U- : 2T‘3T 

2 

ZnGeP 2 

5.49 

10.80 

1.97 

1.8-2.3 (2.10) 

— 

10 13 -10 15 

20 

10 10 -10 17 

10-12 23 27 34 

3 

ZnSnP 2 

5.65 

11.30 

2 

1.6-2.1 (1.86) 

— 

— 

55 

10 16 -10 17 

5|10 12 23 34 [54 

4 

MgSiP 2 

5.72 

10.11 

1.77 

2.2 (2.25) 

— 

— 

— 

— 

26 96 

5 

CdSiP 2 

5.68 

10.44 

1.84 

2.2 (2.10) 

200-1,000 

10 lo -10 15 

— 

— 

11 12;23 26 34 

6 

CdGeP 2 

5.77 

10.82 

1.88 

1.6-1.8 (1.91) 

100 

10 n -10 14 

25 

10 10 -10 15 

10 12 23 31 

7 

CdSnP 2 

5.90 

11.52 

1.95 

1.0-1.5 (1.33) 

2,000 

10 15 -10 18 

— 

10 14 

10 12 23 34] 

8 

ZnSiAs 2 

5.61 

10.88 

1.94 

1.7-2.1 (1.53) 

40 

10 s 

140-170 

10 13 -10 17 

6^ 12 23 26 34 

9 

ZnGeAs 2 

5.67 

11.15 

1.97 

0.6-1.1 (1.21) 

— 

— 

55 

10 16 -10 19 

7 10 12 34 

10 

ZnSnAs 2 

5.85 

11.70 

2.00 

0.6-0.7 (0.89) 

— 

10 15 

300 

10 17 -10 21 

12 23 ; 

11 

MgGeAs 2 

5.66 

— 

— 

- (1-82) 

600 

10 18 

35 

10 19 

1971 

12 

CdSiAs 2 

5.88 

10.88 

1.85 

1.5-1.6 (1.57) 

— 

10 17 

500 

10 14 -10 17 

11 12 23 26.34 

13 

CdGeAsa 

5.94 

11.22 

1.89 

0.5-0.6 (0.70) 

2,500 

10 16 -10 18 

1,500 

10 16 -10 18 

1121341 

14 

CdSnAsa 

6.09 

11.94 

1.96 

0.3 (0.38) 

12,000 

10 17 -10 18 

190 

10 17 -10 18 

12 23 34 98 


ZnSiP 2 Single Crystal XRD Results 

Table S2 Atomic coordinates and site occupancies for ZnSiP 2 . 


Atom 

Wyckoff site 

X 

y 

Z 

Occupancy 

> 

to 

Zn 

4a 

1/2 

0 

3/4 

0.99(2) 

0.005(1) 

P 

8d 

0.73023(6) 

1/4 

1/8 

1.00(2) 

0.005(1) 

Si 

4b 

0 

0 

1/2 

1.00(2) 

0.005(1) 


1-0 I 11 




































Table S3 Crystallographic data for ZnSiPV 


Item [units] 

Value 

Formula 

ZnSiP 2 

Space group 

I42d (No.122) 

Crystal system 

Tetragonal 

a [A] 

5.3986(2) 

c [A] 

10.4502(6) 

v [A 3 ] 

304.57(2) 

z 

4 

FW [g/mol] 

155.413 

Pealed, [g/cm 3 ] 

3.389 

Atomic density [atoms/cm 3 ] 

5.25 x 10 22 

T [K] 

293(2) 

A [A] 

0.71073 

@ maximum 

49.81 

Number unique reflections (n) 

798 

Number reflections I >2cr(I) 

726 

Number of refined parameters (p) 

11 

Extinction coefficient 

0.046(2) 

/i [mm -1 ] 

9.18 

Flack parameter 

0.023(6) 

R(int) [%] 

1.93 

i?(F) a [%] 

1.71 

Rw(F 2 ) b [%] 

3.79 

gof(f 2 ) c 

1.099 

^Pmini ^Pmax 

-1.198, 0.583 


a 

b 

c 


i?(F) = S||F 0 |-|F c ||/S|F 0 | 

Rw(F*) = [Xw(F 2 - F 2 ) 2 /Ew(F 2 ) 2 \ 1/2 

GOF{F 2 )= [(£ w/ \F 2 - F 2 \ 2 )/{n — p) 



Table S4 Anisotropic displacement parameters for ZnSiP 2 - 


Atom 

Uu 

U22 

U33 

Zn 

0.00463(5) 

0.00463(5) 

0.00522(6) 

P 

0.00481(9) 

0.00455(8) 

0.0045(10) 

Si 

0.0046(1) 

0.0046(1) 

0.0042(1) 


Table S5 Selected bond distances (A) and angles (°) for ZnSiP 2 . 


Configuration 

Distance 

Configuration 

Angle 

Zn-P (x4) 

2.3767(2) 

P-Zn-P (x2) 

107.582(3) 



P-Zn-P (x2) 

113.320(7) 

Si-P (x4) 

2.2522(2) 

P-Si-P (x2) 

109.102(7) 



P-Si-P (x2) 

109.656(4) 


In the SiP 4 tetrahedra, P-Si-P angles measure 109.102(7)° and 109.656(4)°, which are close to the ideal 109.5° 
angles expected for tetrahedral geometry. In ZnP 4 , however, two of the angles measure 107.582(3)° while two angles 
measure 113.320(7)°, which indicate a distorted tetrahedral Zn environment. 
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Non-Dimensional Form of DAP Time Dependent Equations 

The non-dimensional variable, x = r 3 N, and parameter, £ = a 3 N, are defined from r, the distance between DAP’s, 
a, the characteristic distance (attributed to the larger of the donor or acceptor effective Bohr radiP®), and A, the 
concentration of the majority defect. The recombination rate as a function of x is then 


W{x) 


Wq exp 


-2 (x/C) 1/3 


(SI) 


where Wq is the rate as r —> 0, as it is in Eqn. [ 2 J Given that x = r 3 N —» dx = 3 r 2 N dr —>■ r 2 dr = the PL 

intensity as a function of time, t, according to Eqn (|3]) , but in terms of x instead of r, is proportional to 


I{t) oc exp 


-4 ^00 

— 7r / {exp[—W{x)t] — \)dx 

y Jo 


4 r °° 

x-n W(x) exp [-W (, x ) t] dx, 
3 Jo 


(S2) 


Intrinsic Point Defect Formation Enthalpies 


Doping of many wide-band gap semiconductors has proven to be challenging.^ Figure [4] (a) shows a heat map of Ep, 
which is determined by the charge balance between the predominant donor and acceptor point defects. The Fermi 
level lies above mid-gap over the majority of the single phase region where ZnSiP 2 is stable, resulting in intrinsic to 
moderately 71 -type conductivity. This result is promising as many compound semiconductors suffer from degenerate 
conductivity. In fact, ZnSiP 2 has been synthesized as both n- or p-type material, likely due to impurities arising from 
the synthesis techniques. Crystals grown in a Zn flux have been reported with both w-type,^ 22 * 25 * 30 * 3 1 ^ and p-type 
conductivity.^ ^ 25 * 31 ! All reports of growth in Sn fl ux ha ve been n-type.& b 22 * 24 * 25 * 29 * 31 * 33 * 36 J Halogen assisted vapor 
transport growth with I has produced n- type crysta ls,while use of Cl as the carrier gas, with ZnCl 2 and PbCl 2 as 
Cl sources, has produced p-type crystals.^ 1 — 25 There is one report of n-type conductivity resulting from Cl vapor 
transport using SiC4 as the Cl source)- 22 Th e crystals ha ve been intentionally, extrinsically doped n-type by adding 
Se, Te, In, or Ga, and p-type by adding Cu. itESEIllIIll] The electronic properties that have been reported in these 


and other studies are given in Table SI When other elements are involved in the synthesis (Sn flux or I, or Cl vapor 
transport) they have been found as impurities in the resultant crystals.^ These studies demonstrate that ZnSiP 2 can 
be synthesized either n- or p-type, using extrinsic dopants, with carrier concentrations acceptable for PV devices. 



E f (eV) E F (eV) E F (eV) 


Fig. SI Defect formation enthalpies (AHd,,) of 7 different point defects (vacancies and antisites) in all possible charge 
states (q ranging from 4- to 4+) at points 1, 2, and 3 on the phase diagram (a). The Fermi level is shown as a vertical dashed 
line along with mid-gap which is shown as a vertical dotted line. 
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Stability of Photoelectrodes 


The stability of n-type semiconductors under aqueous anodic conditions has been a challenge for many photoelectro¬ 
chemical applications (e.g. solar water splitting). The aqueous Fe(CN)g 3 system has been used to monit or the 
water stability of many semiconductor anode materials with and without different protection layers.®- 101 When 
ZnSiP 2 electrodes were tested under aqueous conditions for long periods of time (2-8 hours), some dissolution of 
the semiconductor was observed by scanning electron microscopy (SEM) and the plrotocurrent was observed to de¬ 
crease slightly while the open circuit potential ( E oc ) remained relatively constant. To confirm that the dissolution 
current of the ZnSiP 2 was negligible, an electrode was also tested in a non-aqueous electrochemical cell using the 
ferrocene/ferrocinium (Fc +/, °) redox couple (10 rnM ferrocene, 0.05 mM ferrocinium, LiC 104 supporting electrolyte in 
dry acetonitrile).^ The photocurrent was very similar (within 10%) between the aqueous (50/350 mM Fe(CN)g ) 
and non-aqueous measurements; this is consistent with the expected photon flux through the solution being nearly 
the same (within 1%) in each case. This similarity in photocurrent confirms that dissolution was not a significant 
component in the J-E photoresponse observed for ZnSiP 2 electrodes under aqueous conditions. 
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